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1.  INTRODUCTION 

For  several  years,  the  Weapons  Systems  Research  Laboratory  (WSRL)  has  been 
concerned  with  the  external  carriage  and  separation  of  stores  from  aircraft. 

As  part  of  this  general  field  of  research,  it  is  frequently  necessary  to  carry 
out  wind  tunnel  testing  to  provide  data  on  the  flow  fields  close  to  an  aircraft 
and  so  determine  the  conditions  a  store  is  likely  to  encounter  in  this  region. 

In  addition,  there  is  alwuys  a  need  to  reduce  the  installed  drag  of  external 
stores  because  of  their  considerable  effect  on  aircraft  performance.  This,  in 
turn,  implies  a  need  for  an  accurate  means  of  measuring  the  installed  drag  to 
determine  the  effects  of  any  modifications  to  the  store  or  aircraft. 

This  report  describes  a  facility  which  has  been  developed  at  WSRL  to  carry 
out  these  two  tasks.  The  S-l  wind  tunnel  and  associated  data  processing  equip¬ 
ment  were  employod  together  with  an  IBM  370/3033  computing  system  for  off-line 
data  reduction.  The  facility  consists  of  a  yawmeter  probe  mounted  on  a  computer- 
controlled  traverse  gear  system  in  the  working  section  of  the  wind  tunnel.  Under 
the  control  of  a  PDP-1S  computer,  the  probe  is  moved  through  a  set  of  points  in 
a  rectangular  box  shaped  grid.  The  data  recorded  at  each  grid  point  are  pro¬ 
cessed  through  the  data  acquisition  system  of  the  S-l  complex  and  recorded 
on  magnetic  tape  which  then  undergoes  the  final  phase  of  data  reduction  on  the 
IBM  370/3033  machine  off-line.  This  last  step  enables  CALCOMP  plots  of  desired 
flow  field  parameters  to  be  provided.  Furthermore,  if  the  appropriate  traverse 
grid  has  been  chosen,  the  momentum  defect  in  the  wake  behind  a  model  can  be 
computed  and  used  to  obtain  the  profile  drag. 

The  individual  components  of  the  facility  are  described  in  further  detail  in 
Section  2  of  this  report.  The  off-line  data  reduction  program  is  discussed  in 
Section  3.  Section  4  presents  some  applications  of  this  facility  together  with 
some  experimental  comparisons  aimed  at  validating  the  drag  calculations. 


2.  FACILITY  DLSCR1 PT10N 


2.1  The  wind  tunnel 

The  facility  discussed  here  was  designed  for  use  with  the  S-l  wind  tunnel 
in  the  Acroballistics  Division  of  WSRL.  It  is  a  continuous  flow  tunnel, 
capable  of  subsonic  operation  up  to  Mach  1  and  from  Mach  numbers  of  1.4  to 
2.8  when  used  supersonically.  The  standard  subsonic  working  section  is 
approximately  0.4  m  x  0.4  m  in  cross-section. 

2 . 2  The  probe 

Exploration  of  the  flow  field  is  performed  using  an  hemispherical  headed 
yawmeter  probe  designed  to  measure  Mach  number,  flow  angles  and  dynamic 
pressure.  Geometrical  details  of  the  probe  are  given  in  figure  1,  and 
figure  2  shows  it  in  position  close  to  a  model  in  the  wind  tunnel.  Pressure 
tappings  in  the  probe  are  connected  to  two  differential  and  two  absolute 
pressure  transducers  outside  the  working  section  via  flexible  tubing.  The 
differential  transducers  arc  connected  across  holes  1  to  3  and  2  to  4  (see 
figure  1)  and  so  give  outputs  dependent  on  the  flow  angles  (downwash  and 
sidewash)  experienced  by  the  probe.  The  absolute  transducers  measure  the 
pressures  at  hole  S  and  the  manifold. 

The  electrical  outputs  of  the  transducers  are  fed  via  appropriate  inter¬ 
faces  to  the  data  acquisition  system  attached  to  the  tunnel  and  controlled 
by  the  PDP-1S  computer.  Further  details  on  the  design,  calibration  and 
setting-up  of  the  probe  are  given  in  reference  1,  while  the  software  and 
hardware  concerned  with  the  data  acquisition  system  are  discussed  in  detail 
in  references  2  and  3. 
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2.3  lh<.  tit.vor so  gear 

The  probe  is  attached  to  a  traverse  gear  which  is  capable  of  motion  in  3 
orthogonal  directions  under  either  manual  or  computer  control  by  the  PDP-15 
computer.  Full  details  of  both  the  hardware  and  the  computer  commands 
needed  to  invoke  the  grid  traverse  program  are  also  given  in  references 
2  and  3.  For  the  task  discussed  here,  computer  control  with  the  standard 
grid  traverse  program  was  used  exclusively. 

The  computer  is  programmed  to  move  the  probe  through  a  rectangular  box 
defined  in  terms  of  "box  axes"  by  starting  coordinates  (Xg,Y<.,Zg), 

finishing  coordinates  (Xp,Yp,Zp)  and  increments  (A  Xg,A  Yg,A  .  This 

together  with  the  required  pitch  angle  (0^)  of  the  box  relative  to  wind 

tunnel  axes  (Xj-.Y-pILj.) ,  defines  a  three  dimensional  grid  of  points  with  flow 

field  data  being  measured  at  each  point.  The  axis  systems  and  the  pitch 
angle  0 are  explained  in  figure  3.  Box  axes  and  tunnel  axes  are  assumed 

to  be  based  on  the  same  origin  which,  when  dealing  with  aircraft  models,  is 
frequently  taken  to  be  at  the  centre  of  gravity  of  the  aircraft. 

In  practice,  most  mvestigat ions  employ  only  a  two  dimensional  grid  of 
points.  For  example,  a  wake  traverse  to  compute  the  profile  drag  of  a  model 
requires  a  traverse  grid  in  a  plane  perpendicular  to  the  freestream  (or  Xj.) 

direction.  This  is  obtained  by  setting  Xg  *  Xp,A  ■  0  and  0g  »  0  when 

executing  the  grid  traverse  program  on  the  PDP-1S.  The  majority  of  flow 
field  investigations  have  been  carried  out  in  a  plane  perpendicular  to  the 
Y^.  axis  and  specified  in  an  analogous  fashion. 

Experience  has  shown  that  between  3  and  S  s  settling  time  must  be  allowed 
after  the  probe  is  moved  to  a  new  position  to  enable  the  signals  from  the 
pressure  transducers  to  reach  a  steady  value.  The  grid  traverse  program 
enables  this  delay  to  be  selected  by  the  experimenter  and  it  is  then  auto¬ 
matically  controlled  by  the  computer. 

2.4  The  PDP-15  computer 

As  can  be  inferred  from  the  discussion  so  far,  the  PDP-15  computer  has  a 
major  role  to  play  in  both  controlling  the  execution  of  the  experiment  and 
providing  on-line  data  acquisition.  During  the  running  of  the  experiment, 
real-time  data  processing  is  confined  to  analysing  the  outputs  from  the 
pressure  transducers,  converting  these  to  pressures  using  calibration  data 
held  in  the  computer  and  storing  the  results  on  7-track  magnetic  tape  (the 
primary  tape).  Later,  the  PDP-15  "probe  processor"  program  is  used  off¬ 
line  to  combine  this  information  with  probe  calibration  data  established 
separately (ref . 1)  and  to  store  the  computed  flow  field  parameters  at  each 
grid  point  on  a  second  7-track  tape  (the  secondary  tape).  This  tape  is 
then  submitted  for  processing  on  the  IBM  370/3033  with  program  IG2.61.1/D. 
Subsequent  stages  of  processing  are  discussed  in  Section  3  while  the  general 
flow  of  data  manipulation  is  represented  schematically  in  figure  4. 


3.  PROGRAM  IG2.61.1/D 

3.1  Program  organisation  and  availability 

The  version  of  the  data  reduction  program  IG2.61.1/D  (version  1.4) 
described  here  consists  of  four  subprograms  which  are  executed  by  a  short 
main  program  at  the  discretion  of  the  user.  It  is  currently  coded  in  IBM 
FORTRAN  IV  for  use  on  the  370/3033  computer  operating  in  the  batch  or  non¬ 
interactive  mode.  Each  subprogram  is  logically  independent  of  any  other 
and  can  be  executed  separately.  As  they  are  processed  through  the  sub¬ 
programs,  the  data  originally  generated  in  the  wind  tunnel  are  stored  on 
intermediate  disc  datasets  for  use  by  other  subprograms. 
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These  datasets  are  formatted  to  enable  any  spurious  data  to  be  modified 
manually  by  using  the  IBM  Time  Sharing  Option  (TSO)  EDIT  facility  from  a 
remote  terminal (ref . 4) . 

It  is  hoped  that  most  error  conditions  have  been  anticipated  and  indicated 
by  messages  in  the  printed  output  generated  by  the  program.  A  detailed 
discussion  of  the  data  roquired  to  be  supplied  by  the  user  will  not  be 
given  here  as  the  listings  of  both  the  program  IG2.61.1/D  and  the  Job 
Control  Language  (JCL)  statements  needed  to  execute  it  on  the  IBM  machine 
are  available  from  the  Aerodynamic  Research  Group  in  W’SRL.  The  program 
listing  contains  numerous  "comment"  statements  explaining  the  function  of 
each  data  item  which  must  be  provided.  This  listing  also  includes  comments 
detailing  the  contents  of  each  word  in  each  record  on  the  intermediate  disc 
datasets.  The  program  can  also  be  supplied  in  the  form  of  80  column  punched 
card  decks  or  copies  on  magnetic  tape. 

For  the  convenience  of  the  user,  Table  1  has  been  prepared.  This  relates 
the  symbols  used  in  this  report  to  their  equivalent  variable  names  in  the 
program.  It  is  expected  that  this  list  should  simplify  either  understanding 
the  detaileJ  workings  of  the  program  or  modifying  it. 

In  addition  to  the  main  calling  program  (labelled  MAIN),  1G2.61.1/D  con¬ 
sists  of  the  four  subprograms  (COPY,  DATRED,  PI.OTG  and  DRAG)  mentioned  above 
and  11  other  subroutines.  The  function  of  each  subprogram  will  be  discussed 
in  the  following  sections. 

3.2  Subprogram  COPY 

Each  secondary  tape  produced  by  the  probe  processor  (labelled  logical 
unit  9)  is  first  copied  to  a  disc  dataset  (labelled  logical  unit  10)  and 
this  is  the  task  of  subprogram  COPY.  Data  arc  read  sequentially,  one 
record  at  a  time  off  the  tape.  Each  record  is  checked  to  sec  if  it  is 
within  the  range  of  records  specified  by  the  user  to  be  accepted  and,  if  so, 
copiod  to  unit  10.  Each  record  contains  data  measured  at  one  grid  point. 

No  data  reduction  is  done  at  this  stage  but  experience  has  shown  the  follow¬ 
ing  two  features  to  be  desirable  and  they  have  been  incorporated  in  COPY. 

(a)  Particularly  when  measuring  drag  by  the  wake  traverse  technique,  it 

is  often  realised  after  a  grid  has  been  executed  and  the  data  recorded 
that  it  was  not  large  enough  to  cover  the  required  region  of  the  flow 
field.  Provision  is  therefore  made  in  COPY  to  merge  the  results  from 
an  additional  grid  obtained  from  a  later  tunnel  run  to  data  already  on 
unit  10.  Provided  the  combined  grid  contains  all  the  points  needed 
to  define  a  grid  of  regularly  spaced  points,  subsequent  processing 
will  treat  the  data  now  on  unit  10  as  if  they  were  all  recorded  during 
the  same  run  of  the  grid  traverse  program.  The  order  in  which  the 
data  records  are  copied  to  unit  10  is  not  important. 

(b)  Human  error  can  cause  a  data  word  in  a  record  to  be  recorded  incorrect¬ 
ly  as  the  experiment  is  in  progress.  For  example,  the  experimenter 
manually  sets  a  run  number  and  a  nominal  working  section  Mach  number 

to  be  written  on  each  record  on  the  primary  and  secondary  tapes. 

Those  particular  parameters  do  not  affect  either  the  running  of  the 
experiment  or  any  other  data  recorded  on  tape.  However,  if  they  are 
recorded  incorrectly,  subsequent  processing  might  be  affected.  COPY 
possesses  the  capability  to  change  any  data  word  in  each  record  within 
a  specified  range  of  records  to  any  given  value.  This  can  be  done 
for  any  or  every  record  read  off  the  secondary  tape  and  causes  it  to 
be  corrected  before  it  is  copied  to  unit  10.  "Bulk  editing"  in  this 
fashion  has  saved  a  considerable  number  of  repeated  wind  tunnel  runs. 

3.3  Subprogram  DATRED 

Limitations  of  the  wind  tunnel  data  acquisition  system  mean  that  the 
experimenter  cannot  define  the  spatial  order  in  which  the  grid  is  traversed 
and  the  subsequent  sequence  of  records  on  both  the  primary  and  secondary 
tapes  and  the  dataset  unit  10. 
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This  is  inconvenient  as  later  processing  assumes  that  data  are  available  in 

some  known  order,  e.g.  ,  records  recorded  at  increasing  Xn,  the  YD  and  finally 

b  b 

Zg  grid  coordinates.  DATRiiD  has  the  ability  to  re-arrange  data  read  off 

unit  10  and  store  them  on  another  disc  dataset  (logical  unit  11)  in  any 
specified  order.  To  do  this,  a  parameter  k^  must  be  specified  by  the  user. 

k^  takes  an  integer  value  between  1  and  6  corresponding  to  the  six  different 

ways  of  arranging  the  grid  coordinates  X  YD,Z„.  Table  2  indicates  the 

O  b  b 

value  of  appropriate  to  the  various  combinations.  That  coordinate  allow¬ 
ed  to  vary  most  rapidly  from  record  to  record  on  unit  11  is  in  column  1  of 
the  table  while  the  coordinate  which  changes  least  is  in  column  3.  To 
obtain  the  drag  on  a  body  using  the  wake  traverse  technique,  k^  must  take 

the  value  of  1  or  2.  Most  flow  field  studies  have  been  performed  parallel 
to  the  XgZg  plane  with  k^  equal  to  3  or  4. 

The  sorting  process  depends  on  the  user  repeating  the  information  which 
defined  the  original  traverse  box,  namely  the  starting  coordinates 
(Xs^gtZg),  finishing  coordinates  (Xp.Yp.Zp),  increments  (A  Xg.A  Yg,A  Zfi) 

and  the  pitch  angle  of  the  box  (0g)  which  together  define  a  nominal  sot  of 

grid  points.  Depending  on  the  sequence  indicated  by  k^,  an  attempt  is  made 

to  match  the  coordinates  of  each  point  in  the  nominal  grid  with  those  of  a 
point  on  unit  10.  A  positional  tolerance  (fcv)  is  applied  when  matching  the 
points  to  allow  for  accumulated  truncation  errors  during  the  various  stages 
of  processing.  If  a  match  is  not  made,  an  error  mossage  is  printed  advising 
the  usor  to  cither: 


(a)  repeat  the  whole  tunnol  run  and  all  subsequent  stages  of  data  process¬ 
ing  or, 

(b)  repeat  the  tunnel  test  for  the  missing  point  and  merge  it  to  the 
existing  data  on  unit  10  using  COPY.  hxecute  DATRLD  again,  or 

(c)  insert  a  manually  computed  record  into  the  data  on  unit  10  or  11 
using  the  FDIT  facility  of  IBM  TSO. 


A  complete  set  of  data  at  all  nominal  grid  points  must  be  present  on  unit 
11  for  subsequent  stages  of  processing  to  be  successful. 

During  this  step,  the  measured  flow  angles 


(Qpt  P )  are  used  to  compute 


flow  angles  p^)  in  tunnel  axes 

The  following  expressions  arc  used: 


(XT,YT,ZT  in  figure  3). 


tan 


( 


_  tan  (V _  _\ 

cos(?p)  +  tan(op) sin(flp)  J 


(1) 

(2) 


This  conversion  enables  the  probe  to  be  pitched  by  an  angle 


e 

P 


(see  figure 


3)  to  reach  otherwise  inaccessible  regions  near  a  model,  yet  still  provide 

flow  angle  data  in  a  standard  set  of  axes.  Figure  2  depicts  the  probe 

pitched  towards  the  floor  of  the  working  section  such  that  0  *  -10°. 

v  P 


s 
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3.4  Subprogram  PLOTG 

Flow  field  and  wake  traverse  surveys  invariably  generate  large  quant  ties 
of  data  which  cannot  readily  be  absorbed  if  presented  solely  as  a  listing  of 
data  recorded  at  each  grid  point.  With  this  in  mind,  PLOTG  was  written  to 
take  information  off  unit  11  and  plot  it  on  a  CALCOMP  plotter.  A  typical 

plot  is  given  in  figure  S.  In  this  example,  the  total  pressure  measured 
at  the  probe  divided  by  the  freestream  total  pressure  (labelled  PTOTP/PTOTO 
on  the  plot)  is  plotted  aguinst  the  dimensionless  box  coordinate  Zg/c 

(labelled  IB/CBAK).  c  is  a  non-dimensionalizing  length.  Each  separate 
graph  in  this  particular  plot  pictures  the  results  at  a  constant  value  of 
Yg/c  whose  value  is  given  at  the  top  of  each  graph  (labelled  YB/CBAR). 

Xn/c  is  constant  for  the  whole  plot  and  is  given  in  the  top  left  corner 

(labelled  XB/CBAR) .  The  horizontal  axis  scale  applies  to  the  left  most 
graph  only  and  each  subsequent  graph  is  shifted  to  the  right  by  a  fixed 
amount  designated  by  the  user  who  also  supplies  the  title  printed  across 
the  top  of  the  plot. 

Six  parameters  can  be  plotted  on  the  vertical  axis.  They  are: 

(1)  Downwash,  (O.^.) 

(2)  Sidewash, 


(3)  Local  Mach  number,  (M^) 

(4)  Local  dynamic  pressure/Frecstrcam  dynamic  pressure,  (Hp/M0) 


(5)  Local  total  pressure/Frecstrcam  total  pressure. 


(p  /r  ) 

v  p  oJ 


(6)  Local  static  pressure/} reest ream  static  pressure. 


(yp0) 


The  vertical  axis  represents  one  of  the  following  three  non-dimensional 
positional  coordinates,  (X_/c,  Y'/c,  2n/c) .  The  relevant  one  is  automatic- 

ODD 

ally  selected  by  PLOTG  depending  on  tho  value  of  chosen  when  DATRED  was 

executed  and  it  corresponds  to  the  grid  coordinate  given  in  column  1  of 
Table  2. 

A  plot  of  this  nature  can  only  represent  the  results  of  a  two-dimensional 
traverse,  even  though  the  original  wind  tunnel  traverse  might  have  boen  in 
three  dimensions.  Consequently,  PLOTG  assumes  that  the  data  on  unit  11  is 
divided  into  blocks,  each  block  consisting  of  a  set  of  data  taken  at  a 
constant  value  of  that  grid  coordinate  which  varies  least  rapidly  on  unit  11 
(i.e.,  the  coordinate  in  column  3  of  Table  2).  The  user  must  indicate  the 
number  of  the  block  on  unit  11  to  be  plotted.  One  block  is  plotted  on  each 

call  to  PLOTG  from  the  main  program. 

Each  symbol  on  the  graphs  represents  one  grid  point  and  provision  is  made 
in  PLOTG  to  correct  these  points  either  by  a  straight  line  or  a  smooth  curve 
generated  by  the  cubic  spline  curve  fitting  technique  discussed  in  Appendix  1 


3.S  Subprogram  DRAG 

If  a  two-dimensional  rectangular  grid  perpendicular  to  the  freestream 
direction  downstream  of  a  model  is  explored,  the  resulting  distribution  of 
total  pressure  and  Mach  number  can  be  used  to  estimate  the  profile  drag 
experienced  by  the  model.  This  is  expressed  as  a  drag  coefficient  (C^) 

defined  in  equation  1 1. 3  of  Appendix  II  where  the  analysis  is  presented  in 
some  detail.  The  calculations  are  performed  by  subprogram  DRAG. 
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As  with  PLOTG,  the  data  on  unit  11  are  assumed  to  be  divided  into  blocks, 

each  block  containing  data  from  one  wake  traverse  at  one  value  of  XD.  As 

b 

mentioned  in  Section  3.3,  this  means  that  must  be  specified  as  either  1 

or  2  when  the  data  from  unit  10  is  passed  through  HATRED.  One  block  of 
data  is  processed  on  each  entry  to  DRAG  and  the  user  must  supply  the  number 
of  the  block  on  unit  11  which  is  to  be  analysed. 

4.  LXAMPLES  OF  DRAG  CALCULATIONS 


4.1  General  comments 


To  date,  the  calculation  of  drag  using  1G2.01.1/D  has  been  limited  to 
finding  the  installed  drag  on  various  combinations  of  stores  carried  on 
racks  and/or  pylons  mounted  externally  beneath  the  wing  or  fuselage  of  an 
aircraft.  To  do  this,  a  wake  traverse  of  sufficient  extent  to  include  the 
entire  effects  on  the  wake  of  the  model  attributable  to  the  store  is  perform¬ 
ed.  The  installed  drag  of  the  store  alone  is  obtained  by  repeating  the 
same  wake  traverse  with  the  store  removed  and  differencing  the  two  results. 

The  required  size  of  the  traverse  would  be  difficult  to  ascertain  were  it 
not  for  the  plotting  facility  available  in  subprogram  PLOTG.  A  plot  of 
total  pressure  (as  in  figure  S)  quickly  establishes  whether  or  not  a  grid 
is  large  enough.  It  also  provides  an  indication  of  the  spacing  between  grid 
points  roquirod  to  resolve  variations  in  the  wake  parameters.  A  typical 
wake  traverse  grid  to  measure  the  installed  drag  on  a  store  mounted  on  a 
cylindrical  fuol  tank  beneath  the  wing  of  a  1/24  scale  model  with  a  wingspan 
of  0.343  m  would  be  (in  inches  -  as  required  by  the  grid  traverse  program  on 
the  PDP-15  computor) 


Starting  coordinates:  •  -0.0,  Y<.  »  -0.9,  2,- 


-1.2 


Finishing  coordinates:  Xp  *  -b.0,  Yp  »  *1.0,  2f  ■  *1.2 

Increments:  A  XB  ■  0.0,  A  \  •  *0.1,  A  :  ■  *0.1 

b  B  B 

Box  incidence:  ■  0.0  Jcgrccs. 


The  origin  of  these  coordinates  is  the  rear  of  the  centre  line  of  the  fuel 
tank. 


4.2  Two  experimental  comparisons 

Two  cases  nro  given  which  compare  the  technique  discussed  in  this  report 
to  compute  the  drag  of  external  stores  with  other  experimental  and  published 
values . 

Figure  6  shows  the  experimental  results  from  tests  with  a  1/50  scale  model 
of  a  variable  geometry  aircraft  carrying  six  stores  on  one  rack/pylon  assem¬ 
bly  under  each  wing.  Results  arc  given  for  the  drag  increment  (A  Cp)  due 
to: 

(i)  the  addition  of  the  stores  alone,  and 
(ii)  the  effects  of  carrying  the  bare  rack/pylon. 

These  arc  compared  with  data  derived  from  a  combination  of  flight  and 
wind  tunnel  testing  using  a  geometrically  similar  store  and  published  in  the 
aircraft’s  flight  manual.  Two  cases  of  nominal  wing  sweep  (40°  and  54°) 
were  considered  at  Mach  numbers  between  0.7  and  0.9.  The  differences 
between  the  two  sets  of  results  are  within  the  bounds  normally  accepted  for 
drag  measurements,  especially  after  the  large  differences  in  scale  and 
Reynolds  number  between  the  sources  of  data  are  considered. 
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The  second  example  is  depicted  in  figure  7.  In  this  case,  the  comparison 
is  with  existing  data  obtained  previously  at  WSRL  from  a  wind  tunnel  test  on 
a  1/50  scale  model  of  an  aircraft  carrying  four  stores  on  an  underwing  carrier 
on  each  wing.  A  Mach  number  of  0.8  was  used.  For  the  original  tests,  the 
model  was  mounted  on  a  strain  gauge  equipped  sting  type  force  balance  and  the 
drag  recorded  directly  at  several  aircraft  angles  of  attack  between  -5°  and 
♦10°.  The  wake  traverse  tests  were  performed  on  the  same  model  at  the  same 
Reynolds  number  using  a  simpler  uninstrumented  sting  support .  The  deflec¬ 
tion  of  the  sting  under  aerodynamic  loads  was  taken  into  account  when 
obtaining  the  actual  angle  of  attack  of  the  model. 

There  are  a  number  of  possible  explanations  for  any  differences  between 
the  results  from  the  two  techniques  but  the  most  significant  ones  appear  to 
be  that  the  accuracy  of  the  drag  component  of  the  force  measurement  balance 
is  approximately  *0.0010  in  A  C^,  close  to  the  observed  differences.  In 

addition,  at  the  time  the  force  measurement  tests  were  made,  the  Mach  number 
distribution  in  the  wind  tunnel's  working  section  possessed  several  irregu¬ 
larities,  especially  near  the  model.  This  has  been  since  corrected  by  a 
detailed  re-design  of  the  model  support  system  and  tunnel  walls(ref .5) . 


5.  CONCLUSIONS 

The  wind  tunnel  facility  described  in  this  report  has  reached  a  high  level  of 
operational  efficiency.  It  has  been  used  for  a  number  of  projects  over  a  period 
of  nearly  two  years  and,  wherever  possible,  has  been  verified  by  comparison  with 
data  derived  from  other  sources. 
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N 


P 

S 

xb,yb,1-b 

X  Y  2 

*  r  *  p 

XS,YS,2S 

*pyt,zt 
A  Va  :b 

c 

k 

1 


NOTATION 

reference  area  for  (equation  1 1. 3) 

constants  of  integration  (equations  1.3,  1.4) 
drag  coefficient  (equation  11.3) 

an  increment  in  due  to  the  addition  of  stores 

drag  (equation  11.2  and  figure  8) 
equations  11.7,  11.8  and  11.9 

a  function 
equation  11.5 

an  integral  (equation  1.0) 

an  integral  over  a  sub  interval  (equation  1.8) 

selection  parameter  (see  Table  1) 

the  second  derivative  of  an  interpolation  function 
(equation  1.1)  or  the  number  of  integration  points 
m  the  y  direction,  or  the  Mach  number 

the  number  of  points  in  a  numerical  integration 
in  the  x  direction  (equation  1.7) 

total  pressure 

a  surface  integral 

box  axes  (see  figure  3) 

finishing  coordinates  in  defining  a  traverse  box  (inches) 
probe  axes  (see  figure  3) 

starting  coordinates  in  defining  a  traverse  box  (inches) 
tunnel  .  ^e  figure  3) 

increments  along  box  axes 

the  reference  length  for  coordinates  (inches) 
a  counting  index 

the  distance  between  successive  grid  points 
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P 

4 

u 


Wwp 


W^T 


x.y 


u 


p 


o.°° 

1 

2 


B 

T 

P 

k 


static  pressure 

dynamic  pressure 

velocity  in  the  Xj.  direction 

velocity  components  in  probe  axes  (see  figure  3) 
velocity  components  in  tunnel  axes  (see  figure  3) 


cartesian  coordinates  (usually  with  a  subscript) 
downwash  in  probe  axes  *  tan 


.-1  (  i 


- 1  /  Wt 

downwash  in  tunnel  axes  *  tan  |  — 

Ut 


sidewash  in  probe  axes  •  tan 


1  VJ> 


sidewash  in  tunnel  axes  *  tan 


-1  VT 


ratio  of  specific  heats  ■  1.4 

the  spatial  tolerance  of  a  grid  point 

the  pitch  of  the  traverse  box  (sec  figure  3) 

the  pitch  of  the  probe  (sec  figure  3) 


air  density 

subscripts  (where  not  stated  explicitly  elsewhere) 
pertaining  to  frecstream  conditions 
a  quantity  at  station  1  (figure  8) 
a  quantity  at  station  2  (figure  8) 
a  quantity  at  station  »  (figure  8) 
pertaining  to  the  traverse  box 
pertaining  to  the  tunnel 

pertaining  to  quantities  measured  by  the  probe 
a  counting  index 
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TABLE  1.  SYMBOLS  WITH  EQUIVALENT  VARIABLE  NAMES 


Synbol  in  text 

FORTRAN  variable  nane 

A 

REFA 

C 

CONSTC 

CD 

CD 

K 

CONSTK 

S 

KODORD 

0 

PDRAC 

Fi.f2.f3 

FI ,F2,F3 

G1 ,G2,G3 

RR1 ,  UU1,  UJUI 

I 

SPLINT 

M 

MN,  MACHI 

VM2 

RMACH,  MACH 

:  po*  p~ 

PTOTO 

p  p 

P*  2 

PTOTP 

XB’^  B’^B 

TRX,  TRY,  TRZ 

XEND,  YEND,  ZEND 

XSTART ,  Y START,  ZSTART 

A  VA  h 

XINC,  YINC,  ZINC 

c 

CBAR 

Pq'Poo 

STATO 

|  PP’P2 

STATP 

^O'^00 

QO 

^p’^2 

Q 

1 

a 

P 

ALFA,  DOWN 

°T 

ALFW 

'p 

BETA,  SIDE 

^T 

BETAW 

y 

GAFMA 

E 

ERROR 

% 

THETAB 

e 

P 

THETAP 
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TABLE  2.  THE  VALUE  OF  SEQUENCING  PARAMETER  ^ 


1 

MOST  RAPIDLY 

LEAST  RAPIDLY 

*» 

VARYING  COORDINATE 

COLUMN  2 

VARYING  COORDINATE 

COLUMN  1 

COLUFW  3 

Figure  1.  General  arrangement  of  yamaeter  probe  and  mounting 


-  15  - 
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Figure  3 


Figure  3.  Definition  of  axes  and  pitch  angles 
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Figure  4 


With  The 

rnr-is 

* 


Figure  4.  Flow  of  data  through  the  facility 


Figure  5.  Typical  CALCOMP  plot  from  PLOTG 
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Figure  6 
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Figure  6.  Drag  increments  of  external  stores  -  comparison  of 
experimental  and  handbook  values 
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Figure  7 


Figure  7.  Drag  increments  of  external  stores  -  comparison  of 
wake  traverse  and  force  measurement  techniques 
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APPENDIX  I 

THE  CUBIC  SPLINE  CURVE  FIT  AND  INTEGRATION 


1.1  Ceneral 

The  computer  program  1G2.01 . 1/D  described  in  the  body  of  this  report 
requires  an  algorithm  which  will  provido  an  analytical  expression  for  an 
accurate  curve  through  a  given  set  of  points.  These  curves  are  used  in 
two  roles: 

(1)  to  provide  interpolation  points  between  measured  grid  points  as 
required  by  the  plotting  subprogram  PLOTG.  Because  the  CALOMP 
plotter  is  constrained  to  draw  a  straight  line  between  two  points, 
the  improssion  of  a  smooth  curve  is  generated  by  locating  a  number 
(typically  10  to  20)  of  interpolation  points  between  each  plotted 
grid  point  and  joining  these  with  straight  lines,  and 

(2)  to  form  the  basis  for  a  numerical  surface  integration  over  a 
rectangular  region  as  used  to  evaluate  the  drag  on  a  body  using  the 
wake  traverse  technique.  Equation  1 1. 2  of  Appendix  II  gives  the 
expression  to  be  evaluated  by  this  means. 

Both  rolos  were  filled  by  a  cubic  spline  curve  fit  as  described  and 
dcrivod  by  Walsh  ct  al.(ref.b).  The  following  is  an  extended  version  of 
that  derivation  with  emphasis  on  thoso  specific  aspects  used  in  IG2.61.1/D. 

1.2  The  interpolation  function 

Suppose  a  function  (F(x))  is  prescribed  by  y^  at  N*1  discrete  values  of 
x.  donoted  by  x^  where  k*0,  1,2.  .  N.  Let  denote  the  second  derivative 
of  the  approximate  interj>olat ion  function  y(x).  Further,  assume  that 
varies  linearly  between  any  two  points  so  that 


y  (*) 


*  “k 


x  -x 


k-1 


1 


k 


Vi* 


X  <  x. 


1.1 


l 


where 


Integrating  1.1  twice  and  applying  the  conditions  that 


x  -  x. 


at  y  (xk)  -  yk 


and 


-  v,  at  y(xk-i}  "  yk-i  • 


leads  to 


/(*)  •  — 5 J 


Vl  (*k-*>*  ^  (*-»k-l>f 


(>(  i 


♦  Cx  ♦  K  , 


1.2 
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where 


i*k  -  Xk-i)  «k  -  \-0 


fyk  \  lk)  (hzi  Vi_V) 

"  *xk-i  \t  '  “nrv  *  xk  \  <  *  6  )' 


Equation  1.2  is  the  interpolation  function  used  in  tho  first  role  mentioned 
above.  Tho  only  terms  not  given  explicitly  by  the  set  of  points  (x^ *>,|t )  in 

equations  1.2,  1.3  and  1.4  are  the  M^,  terms.  These  are  obtained  by 

applying  tho  condition  of  continuity  of  slope  at  each  point,  that  is 

y*  (x‘)  -  yx  (x*)  .  1,5 

This  condition  leads  to  a  set  of  linear  simultaneous  equations  which  are 
readily  solved  by  matrix  inversion. 

1.3  Tho  surface  integral 

If  it  is  assumed  that  a  one  dimensional  integral  of  the  type 


F (x)dx 


can  be  represented  by  sums  of  individual  integrals  ( I over  each  sub-interval 


such  that 


k  «  n  k  -  N 

■  •  1  ■■  •  i.  /:::: 


F(x)dx  , 


k  ■  1 


k-1  x  “  xk-l 


it  is  clear  that  the  integral  (equation  1.6)  can  be  evaluated  by  examing  the 


1^  terms  in  isolation. 


fx  ’  *k 

>  •  J,.„. 


F(x)dx 


can  be  approximated  by  the  interpolation  function  derived  above, 
(equation  1.2)  such  that 


x  ■  x. 

q-  / 

J  X  ■  X. 


y(x)dx  . 
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Substituting  equation  1.2  yields 

'k“  JJ  *V  •  Vi)  •  f  («k  *  *k-i)  'k  *  *'k  •  *■“> 

which  can  be  evaluated  since  all  terms  are  known  or  can  be  derived  using 
methods  discussed  in  Section  1.2.  The  approximation  to  the  integral  I 
(equation  1.6)  follows  from  equation  1.7. 

A  numerical  surface  integration  over  a  plane,  rectangular  surface  defined 
in  cartesian  coordinates  (x,y)  and  containing  a  grid  of  points  disposed 
rectilinearly  in  the  x  and  y  directions  merely  involves  the  same  procedure 
carried  out  in  two  directions. 

Consider  the  surface  integral  S  given  by 


S 


*N  ) 

F(x)  dx  (  dy 


1. 11 


Numerical  approximation  of  S  reduces  to  evaluating  the  "inner"  integral 
(within  the  brackets  {  (  )  along  all  the  rows  points  at  each  value  of  x  and 
storing  the  results.  These  form  the  basis  for  a  new  set  of  data  which  are 
then  integrated  by  the  same  technique  in  the  y  direction. 
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APPENDIX  II 

DERIVATION  AND  EVALUATION  OF  THE  DRAG  EQUATION 

The  following  derivation  is  an  extension  of  the  work  of  Schlichting(ref . 7, 
chapter  25).  Lock(ref.8)  can  also  be  consulted  for  a  similar  analysis. 

Consider  three  planes,  all  of  which  are  aligned  perpendicular  to  the  free- 
stream  direction.  These  are  indicated  in  two  dimensions  in  figure  8  and  are 
designated  stations  °°,  1  and  2.  The  first  plane  (station  °°)  is  sufficiently 
far  upstream  of  the  body  whose  drag  is  being  measured  to  be  regarded  as  having 
froestream  conditions  existing  at  all  points  on  it.  Station  1  is  located  far 
enough  downstream  of  the  body  that  the  static  pressure  (pi)  there  is  essentially 
that  of  the  freestream,  i.e.. 


P>  ■  P»  •  11.1 

Typically,  station  1  is  so  far  downstream  of  the  body  that  wake  traversing  is 
not  feasible  there.  A  third  plane  (station  2)  is  therefore  used  for  the  actual 
traverse.  The  static  pressure  at  a  point  on  station  2  (p3 )  is,  in  general 
different  from  p^,  . 

Momentum  considerat ion; ,  and  the  assumption  that  the  total  pressure  along  a 
streamline  between  stations  1  and  2  is  constant  (i.e.,  P2  =  Pj )  leads  to  the 
following  expression  for  the  profile  drag  (D)  on  the  body 


D 


(U**,  -  uf  )  dSi 


11.2 


and  ui  are  the  air  density  and  velocity  in  the  freestream  (X)  direction  at 
station  2.  u^,  and  u,  arc  the  X  velocity  components  in  the  freestream  and  at 

station  1  respectively.  In  terms  of  a  drag  coefficient  (Cjj) 


C 


D 


’}P 


Ubo 


II. 3 


1  [[  fii 

*11  Ko 


11.4 


where  A  is  a  convenient  reference  area.  The  output  from  the  probe  processor 
program  on  the  POP- 15  consists  of  the  freestream  values  of  total  pressure  and 
Mach  number  (P^  and  with  the  same  quantities  measured  at  grid  points  at  the 
measurement  station  2.  (Pj  and  Mi).  These  terms  must  be  used  when  evaluating 
equation  1 1. 4. 

It  is  convenient  to  define 


G, 


Gi 


vm 

Ut¬ 


il  .5 


Using  several  of  the  standard  relationships  of  compressible  gas  dynamics  and  the 
assumptions  that  P}  ■  P, ,  pi  ■  p^  it  follows  that 


■( 


h 

Pi 


11.6 
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Hence  define 


F, 


(7-1) 


II. 7 


in  which  7  is  the  ratio  of  specific  heats  (■  1.40).  In  addition,  define 


Fj 


II. 8 


F, 


II. 9 


It  is  well  established  in  the  literature  that 


\1  ♦  %”f\ 


i 

7-1 


II.  10 


and 


II.  11 


The  lack  of  a  subscript  of  p  u  and  M  implies  a  quantity  measured  at  any  point 
in  the  flow  field.  Referring  specifically  to  stations  1  and  2,  equations  II. 10 
and  II. 11  can  be  expressed  as 


G,  - 

1 

•  1 

1 

7-1 

11.12 

L  1 

• 

» 

9 

Gj  ■ 

Fa 

♦  7-1 

11.13 

Fj 

♦  7-  l_j 

» 

% 

Gj  • 

Fi 

♦7-1.1 

_  F, 

♦  7-1 J 

• 

11.14 

Equation  II. 4  can  now  be  written 

as 

s  • 

fj  G,  G,  (1  -Gj  )  dS,  . 

II. IS 

The  integrand  of  this  equation  is  readily  evaluated  at  each  grid  point  in  the 
traverse  and  the  numerical  integration  technique  discussed  in  Section  1.3  of 
Appendix  I  can  then  be  used  to  compute  Cp. 
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